Conductance scaling at the band center of wide wires with pure non— diagonal disorder 
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Kubo formula is used to get the scaling behavior of the static conductance distribution of wide 
wires showing pure non-diagonal disorder. Following recent works that point to unusual phenomena 
in some circumstances, scaling at the band center of wires of odd widths has been numerically 
investigated. While the mean conductance shows a decrease that is only proportional to the inverse 
square root of the wire length, the median of the distribution exponentially decreases as a function 
of the square root of the length. Actually, the whole distribution decays as the inverse square root 
of the length except close to G — where the distribution accumulates the weight lost at larger 
conductances. It accurately follows the theoretical prediction once the free parameter is correctly 
fitted. Moreover, when the number of channels equals the wire length but contacts are kept finite, the 
conductance distribution is still described by the previous model. It is shown that the common origin 
of this behavior is a simple Gaussian statistics followed by the logarithm of the E = wavefunction 
weight ratio of a system showing chiral symmetry. A finite value of the two-dimensional mean 
conductance is obtained in the infinite size limit. Both conductance and the wavefunction statistics 
distributions are given in this limit. This results are consistent with the critical character of the 
E = wavefunction predicted in the literature. 

PACS numbers: 72.15.Rn, 71.55.Jv, 11.30.Rd 
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I. INTRODUCTION 

The existence of peculiar properties that differentiate 
pure non-diagonal disorder models from disorder models 
including diagonal disorder has a long history. Proba- 
bly, the first contribution along this line corresponds to 
Dyson's work an a one-dimensional phonon model pub- 
lished in 19531 The existence of a divergent density 
of states at the band center of such disordered systems 
implying a -divergent localization length was pointed in 
later workaa'u. Since then many works on this subject 
have been published (a representative list can be found 
in Ref.([|)). The existence of some kind of derealization 
transition at the band center of such models certainly 
disturbs the widely accepted statement saying that the 
specific form of disorder does not matter in single param- 
eter scaling theorycl. 

A bipartite lattice is a lattice that can be divided into 
two sublattices such that the Hamiltonian changes sign 
under a transformation that changes the sign of the wave- 
function on one sublattice. When pure non-diagonal dis- 
order is considered on a bipartite lattice, the electron- 
hole symmetry of the spectrum is not destroyed by dis- 
order. This property has important consequences as, for 
example, the existence of an eigenstate at E = for any 
disorder realization of a system constituted by an odd 
number of sites (the spectrum shows ±E a pairs plus a 
state at E = 0). Recently, several works pointing to 
the exotic behavior of transport properties of quantum 
wires showing chiral symmetry have been publishedBtl 
For example, scaling of the conductance strongly depends 
on the parity of the number of channels along the wire. 
Also, related activity in field theory has produced several 
models demonstrating a derealization transition in the 



vicinity of the zero energy stated. 

In this paper, the scaling properties of the simplest dis- 
order model preserving chiral symmetry have been care- 
fully analyzed at the band center of quantum wires. Nu- 
merical simulation has been used to get the whole dis- 
tribution of conductances. While the mean value of the 
conductance decays algebraically for wires of odd width, 
an alternative measure of the central value of distribu- 
tions -the median- shows an exponential decay. Actually, 
the larger part of the conductance measurements has an 
exponentially small value for large wires. Therefore, al- 
though the conductance distribution is certainly peculiar 
at the band center, I would still use the term exponen- 
tial localization when referring to the scaling behavior 
of the conductance of long wires. Results change when 
width and length of the wire coincide. In this case, the 
numerical simulation presented in this paper uses a dot 
geometry point of view keeping the size of the contacts 
finite while dot area scaling proceeds. The conductance 
distribution converges to a well defined limit that is com- 
patible with the predicted critical behavior of the state at 
the band center. It is shown that the analytic form used 
to fit conductance distributions comes from a new under- 
lying wavefunction statistics describing the distribution 
of weight ratios of the E = wavefunction. 

The format of the paper is as follows. Section II defines 
a quite simple chain model that allows some analytical 
results and an unbound numerical simulation. Section 
III gives a more general disorder model on finite rectan- 
gular clusters of the square lattice. The way in which 
the conductance is calculated is presented in Section IV. 
Numerical results are given in Section V, first for wide 
disordered wires and, second, for square clusters. The 
last Section of the paper compiles the main conclusions 



reached by this numerical study of conductance scaling. 



II. TOY MODEL 

Let me begin with a detailed description of the scaling 
properties of the conductance of a chain showing pure 
non-diagonal disorder. In this case, sign changes of the 
hopping parameters do not matter and changes in their 
absolute value should be considered. One case that al- 
lows some analytic results consists of a chain with hop- 
ping parameters that randomly take values of 1 and V2 
with equal probability. Starting from the selfenergy of an 
ideal chain at its band center E = — i, successive self- 
energies can be obtained through the disordered part of 
the chain by means of the usual iterative sequence (see, 
for example, Eqs. (7)-(10) in Rcf.(|)): 
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where £„ 2 takes values 1 or 2 with equal probability. Since 
selfenergy remains purely imaginary, a more convenient 
form is possible: 



-•n+l 



/i 



(-£»/*) 



(2) 



The conductance G of a sample is obtained from the self- 
energy at the end of the disordered chain part E^r/i: 
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(3) 

Eq. (|3|) shows that the conductance varies between 
and 1 as it corresponds to a one-channel system. Actu- 
ally, the conductance takes the value 1 only if — Ejv/i = 
— Eo/j = 1. Repeated use of Eq.(@) shows that the form 
in which the random hopping elements appear in the self- 
energy expression is: 
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(4) 



Therefore, a perfect transmission through the chain is ob- 
tained when the number of normal (hopping equal to 1) 
and strong (hopping equal to v2) bonds of the numerator 
coincides with the corresponding numbers of the denom- 
inator. Solving this simple combinatorial problem, one 
gets a probability 
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for the peak of .the probability distribution of the conduc- 
tance at G = 10. Notice that the distribution is a sum of 



delta functions since the hopping takes just two different 
values; for example, the peak below the one at G = 1 
appears at G = 8/9 and corresponds to — Eat/z = 2 or 
— Ejv/i = 1/2 (one extra strong bond either in the nu- 
merator or the denominator). An asymptotic expansion 
for large enough chain lengths can be found for the sum 
in (§: 



p(G=l) 



(6) 



The scaling behavior of this peak is enough to explain 
an inverse square root law for the mean value of G given 
that other peaks of the distribution decay in a similar or 
faster way with the size of the chain. This is an inter- 
esting behavior since it perfectly coincides with scaling 
predictions for wide wires of odd number of transversal 
modestl. 

Numerical simulations can be done as precise as nec- 
essary for this simple model. This fact allows a detailed 
comparison with theoretical predictions. Fig. [fl shows 
the scaling of the conductance as described by its mean 
value. Root mean square deviation is about some tenths 
of the conductance unit. In any case, the average con- 
ductance at large sizes (N > 1000) shows a power law 
scaling with the distance that is compatible with the in- 
verse square root law given by Eq.(|6|). According with 
this result, a disordered chain with pure non-diagonal 
disorder shows non-standard scaling at the band center. 
Nevertheless, alternative measures of the central value of 
the distribution restore to some way the usual exponen- 
tial decay of one-dimensional conductances, even in the 
presence of chiral symmetry. For example, the geometric 
meanEj of the same conductance distributions shows a 
much more pronounced decrease with length (see Fig. 0) 
although the corresponding standard deviation is of the 
same order as the mean making doubtful its statistical 
relevance. But there are other alternatives which do give 
a good description of the overall scaling of the distribu- 
tion. Both the medianEiJ or any definition based on the 
value of the integral of the distribution between and an 
arbitrary upper limit G max flow to exponentially small 
values. The physical meaning is clear in this case: half 
or more of the measures are exponentially small at large 
chain lengths. Actually, the precise scaling law for these 
central value alternatives is: 



Gq ~ exp 



(7) 



where £ gives a measure of the exponential localization 
length. Fig. [2] shows that fits according to this law are 
excellent over the whole length range. The ultimate rea- 
son for such statistically disappointing results is simply 
the unusual size scaling of the distribution. While the 
major part of the distribution below G = 1 decays in an 
algebraic form, the weight of the distribution accumu- 
lates in an exponentially small region near G = (Fig. 



[| in Section V illustrates graphically this behavior) . In 
this way, the upper part of the distribution dominates 
the scaling behavior of mean, root mean square devia- 
tion, etc., while accumulation at G = gives the be- 
havior of central value definitions based on the integral 
of the distribution. In my opinion, these last definitions 
are better suited for characterizing the whole distribution 
than standard averages. Ultimately, one should look at 
the precise experimental protocol followed to get a value 
of the conductance before making predictions about the 
result of the measurements. 



take place owing to the existence of a true eigenstate at 
exactly E — for any piece of the system showing an odd 
number of sites (the matrix E — H is locally singular). 
Nevertheless, numerical calculation proceeds straightfor- 
ward when pivoting over the whole Hamiltonian matrix 
that includes the ideal leads is allowed. 



V. RESULTS 



A. Wide disordered wires 



III. NON-DIAGONAL DISORDER MODEL 

The lattice Hamiltonian describing random hopping on 
a L x M cluster of the square lattice is : 
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where c\ creates an electron on site I, I and I' are nearest- 
neighbor sites, and tw is the hopping energy from site I to 
V . It takes values 1 and -1 with equal probability. Let me 
refer to this model as the Random Hopping Sign (RHS) 
model. Obviously, the square lattice can be divided into 
two sublattices such that atoms belonging to one of them 
hop only to sites belonging to the other sublattice when 
described by Hamiltonian (||). Therefore, this Hamil- 
tonian changes sign under a transformation that changes 
the sign of the electron operators on one sublattice. Con- 
sequently, the spectrum is symmetric relative to the band 
center at E — for any disorder realization, i.e., for any 
values of the random variables of the model {tu>}. 

The model given by Eq.(|8J) is probably the simplest 
two dimensional model showing chiral symmetry. Many 
changes can be done to this model preserving chirality. 
For example, absolute values of the hopping can fluctu- 
ate in addition to their signs or complex values of the 
hopping parameters can be considered if random mag- 
netic fluxes are simulated. The results described below 
are not sensible to any of these modifications of the disor- 
der model. Numerical values are somewhat changed but 
trends remain exactly the same. 



IV. CONDUCTANCE CALCULATION 

Randomly generated samples of L x M clusters are con- 
nected to ideal leads of width M. Typically, L > M in a 
wire geometry. The conductance of the whole system is 
obtained using Kubo formulaEEl within exact one-electron 
linear response theory. Computational details are given 
elsewhereEjilj. Let me just mention that the inversion of 
the Hamiltonian matrix needed to get the Green function 
of the wire cannot proceed slab by slab as it is usually 
done within an optimized code. Numerical divergences 



The first aim has been the recovery of some important 
results obtained for wide wires by Mudry, Brouwer, and 
Furusakicl. In particular, the exotic dependence of the 
scaling law on the parity of the wire width is obtained 
for the present model (a simplification of their Random 
Flux Model) . Open boundary conditions have been used 
to get conductances of stripes of fixed width and number 
of open channels (the wire width is equal to the number of 
channels at the band center). Sample lengths have been 
varied from 99 to 1980 in steps of 99. As many as 10 4 
samples are necessary to get good values of means and 
other central values of conductance distributions. Fig. 
shows the scaling law for two typical odd widths (9 and 
19 channels) in a log- log plot. Error bars are comparable 
to the symbols representing the conductance averages. A 
power law fit to the numerical data is compatible with 
a mean conductance scaling proportional to the inverse 
square root of the sample length L: 

<G>r 



L 



This is precisely the scaling law obtained in Ref.ra) for 
quantum wires of an odd section. Notice that scaling pro- 
ceeds smoothly without distinguishing odd and even wire 
lengths. As noted by these authors, the variance of the 
conductance distribution is as large as its mean, making 
the mean a poor characteristic of the whole distribution. 
Actually, the predicted relationship 



<G 2 > 
<G> 



is accurately reproduced by the numerical data (see filled 
symbols of Fig. [| and the dashed lines that are simply 
2/3 of the previous fits). On the other hand, when an 
even number of channels is studied, exponential scaling of 
the average conductance together with an exponentially 
small typical deviation is obtained in complete agreement 
with previous workEl. 

A deeper understanding of scaling of the conductance 
mean value can be obtained from the analysis of the 
whole conductance distribution. Here, conductances for 
10 5 samples of width 9 and lengths 10 2 , 10 3 , and 10 4 have 
been compiled and the corresponding histograms plotted 
in Fig. [|. It is clear that the probability for large values 
of the conductance diminishes as the length of the wire 



increases. Actually, since the plot is semilogarithmic the 
roughly equal separation between filled and empty circles 
for one side and between empty circles and diamonds for 
the other side, implies a power law decrease of the proba- 
bility. The weight that the probability distribution loses 
for large conductance values goes close to G = 0. I have 
checked that the divergence of the probability at G = 1 
is proportional to l/\/l — G while the divergence at the 
origin looks also algebraic but with an exponent starting 
close to — i for small disorder and decreasing towards — 1 
as disorder increases. 

Although the conductance corresponding to ideal (non 
disordered) wires equals their widths (9 and 19, in this 
case), Fig. || shows that conductance just reaches the 
value of 1 for disordered systems. It seems that just one 
channel is effective. Therefore, it is tempting to com- 
pare these conductance distributions with the one cor- 
responding to one random channel within the orthogo- 
nal universality classed arguing that before localization 
the effect of disorder is just randomizing transport coef- 
ficients. But the comparison is very bad since apart for 
small differences all theories show probabilities that con- 
tinuously decrease from G = to G = 1. In particular, 
the simplest result that applies to a random channel de- 
scribed by scattering matrices of the Circular Orthogonal 
Ensemble is: 



C 



I M 2 + M-2 
AL M-l 



(11) 



with I the mean free path, L the wire length, and M the 
wire width (which coincides with the number of channels 
at the bandcenter). A second scaling parameter -q that 
characterizes the disorder on a microscopic scale does not 
appear because it vanishes for the random flux model 
(RFM) and the model studied hep. (see Section III) is 
just a special case of RFM modelliil. When C is small 
enough (C < 0.01) the distribution is very well described 
by a much simpler expression (error smaller than one 
percent for G ^ 0): 



p(G) 
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(12) 



which allows the evaluation of the mean and the variance 
of the distribution: 



<G>~2a/- 
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p(G) = 



2VG 



(9) 



that is quite different from the one obtained by numeri- 
cal simulation. There is a direct mathematical reason ex- 
plaining this failure. Typically, several states contribute 
to the Green function calculated at an arbitrary energy 
within the spectrum of a disordered system. But this is 
not the case when the Green function of a chiral system 
is calculated at E = which is an eigenenergy of the iso- 
lated system with an odd number of sites. In this situa- 
tion, both the Green function of the isolated finite system 
and the one corresponding to the extended system includ- 
ing the leads (and related to the previous one by a Dyson 
equation) are dominated by the pole at E — 0, that is, 
are completely determined by the E = eigenfunction. 
In the next subsection, I will exploit this feature to ana- 
lyze the conductance distribution as a consequence of a 
precise wavefunction statistics. Meanwhile, let us return 
to the theory of Ref. (ph to analyze the numerical results. 
Mudry, Brouwer, and Furusaki give the following ex- 
pression for the conductance distribution of wires with 
an odd number of channels in the localized regime (see 
Eq. (4.10) of the second paper of Ref. (§)): 



p(G) 
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(10) 



where 



< G 2 >= - < G > 

Using the explicit form of C (Eq.rtnh) an expression is 
obtained that can be used to fit the remaining parameter 
I once an enough number of widths and/or lengths have 
been studied: 



<G> 2 L~ — 

IT 



M 2 + M - 2 
M-l 



(13) 



Fig. (ph shows a fit to the mean and variance values 
obtained for sets of 10 4 randomly generated samples of 
length L = 2000 and width M from three to 37. While 
Eq.(|l3|) works reasonably good for large widths it fails 
near M — 3. Actually, an alternative fit by a linear law 
works sensibly better: 



< G > 2 L = aM + b 



(14) 



Consequently, the distribution given by Eq. ( |To[) will be 
used with 



C 



0.3515A/ + 2.3050 



(15) 



to describe the results obtained by numeric simulation. 

Probability p(G) is shown in Fig. Has continuous lines. 
Although the three numerical distributions are nicely re- 
produced, accuracy is better for wider wires as expected 
from the theory (M <C L is assumed in the theory leading 
to Eq. (|l0|)). The divergence at G — 1 is of the inverse 



square root form while the apparent G _1 non-integrable 
divergence at the origin is regularized by the complex nu- 
merator which vanishes at G = 0. In conclusion, Eq.(J10|) 
is a very good description of the numerical data once the 
constant C is properly estimated. 

Let us try an alternative way of characterizing the cen- 
tral value of distributions of the form shown in Fig. y. 
Previous experience with disorder models of this kind 
proves that the median is more robust than the mean for 
some distributions showing large G variances but smaller 
ln(G) variances, i.e., when geometric mean does it better 
that the usual arithmetic meanly. Unfortunately, the dis- 
tribution of ln(G) is also very broad in the present case. 
Scaling results for the toy model of Section II suggest 
the use of the median. Fig. g] gives the median scaling 
obtained for the previously collected conductance distri- 
butions. In can be seen that median values are exponen- 
tially smaller than mean values at large wire lengths. The 
fact that the median scales towards can be inferred from 
the scaling of the histograms given in Fig. [| although, 
the accumulation near exponentially small values of G is 
not visible in the Figure. Even the scaling law of medi- 
ans is the same that describes the chain model (see, Eq. 
(0)). Alternatively, the exponential decay of the median 
as a function of the square root of the wire length can be 
inferred from Eq. (O) once the distribution is written 
as a logarithmic normal distribution using the variable 
change described in the next subsection. 

The practical meaning of the result is quite clear, more 
than half of the conductance observations are exponen- 
tially decreasing as the wire length increases. Actually, 
the fraction of the samples showing exponentially small 
conductances increases with wire length because proba- 
bility below G — 1 decreases monotonically. Just the dif- 
ference with a conventional one-dimensional scaling (as, 
for example, the one obtained for wires of even width) 
comes from the power law decrease of the upper part 
of the probability distribution that should be compared 
with the exponential decrease characteristic of standard 
scaling. 



cavities. This limit is described in a first approximation 
by scattering matrices of the Circular Orthogonal Ensem- 
ble (COE)t2l. Roughly speaking, a conductance about \ 
per channel can be expected. Second, Anderson exponen- 
tial localization would imply an exponentially small value 
of conductances for large enough dot sizes. This limit 
applies to diagonal disorder, for example. At the band 
center of a system with chiral symmetry, numerical sim- 
ulation shows a behavior similar to the ballistic one (see 
Fig. 0). Mean conductance converges to a well-defined 
finite limit while the whole conductance distribution is 
perfectly described by Eq. (JO]) . Since now the system is 
not quasi-oncdimcnsional as wide wires are, one is forced 
to conclude that there should be deep general reasons for 
the validity of Eq. (|l(j) in this context. 

Let me briefly describe the numerical procedure. The 
transmission between any couple of points within the dot 
is calculated and the corresponding conductance distri- 
bution obtained. To this end, two clean infinite chains 
(the leads) are attached through two arbitrarily chosen 
lattice sites within the disordered square sample (the 
dot). For this geometry, the transmission from site r 
to site r' is given by the following expression: 



T 



4 V(r) 2 ip{r') 2 
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(16) 



where ip(j) is the wavefunction at the band center. From 
this equation, < T < 1, and the conductance is be- 
tween zero and one quantum unit for this numerical sim- 
ulation. While the large number of transmission evalua- 
tions (~ \L A ) allows for a very precise calculation of the 
whole conductance distribution, the dependence of the 
conductance on the separation between point contacts is 
not given by the procedure. The eigenstate at the band 
center is obtained by direct inversion of the Schrodinger 
equation (Hip = 0) where H is given by (0). Once the 
E = wavefunction statistics is known, Eq. (|lj) can be 
used to get a conductance distribution. For example, if 
the Porter-Thomas formta were valid: 



B. Two dimensional system 



/(*) 



: exp(-t/2) 



(17) 



While a standard study of the scaling of the conduc- 
tance in 2D would imply the calculation of conductances 
of increasing L x L samples connected to ideal leads of 
width equal to the square side L, in this work I have 
used a dot setup specially designed for the study of just 
one conducting statet-3. Certainly, when the number of 
incoming channels is fixed by a point contact geometry, 
the only factor that affects the value of the conductance 
is the size (area) of the dot. In this way, the presumably 
increase of the conductance due to wider contacts does 
not obscure the underlying scaling law strictly due to the 
increased size of the system. Two different limits are well 
known for large values of L. First, ballistic transport 
through the sample can occur as it happens in chaotic 



where t = Nip(r) being N the number of sites, the con- 
ductance distribution would be given by: 
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which can be integrated to give the final result: 



Ppt(G) 



* y/G(l - G) 



(18) 



(19) 



Although probability ppx(G) reproduces some features 
of the conductance distributions shown in Figs. and R 



for the smaller cluster sizes (for example, the square root 
divergence at G = 1), it is clearly non comparable to the 
accurate result given by Eq.(lO). What comes as a sur- 
prise is the fact that a logarithmic normal distribution 
of the wavefunction ratios t/t' exactly gives the conduc- 
tance distribution proposed by Mudry et aln. That is, 



assuming 



I \ / ^ , G 2\ 
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(20) 



being x — hx(t/t'), the conductance distribution is given 
by 



*»9{*)S{G- * ) 



(21) 



&- + e~" + 2/ 
This integral is easily solved giving p(G) of Eq.(|lO|) with 



G = 



or equivalently, 



arcosh(G" 1/2 ) 



This result gives some clue over the complicate G depen- 
dence that happens to appear in its distribution func- 
tion. Numerical simulation (see Fig. H) shows that 
Eq.(p0|) accurately describes the wavefunction squared 
ratios of large disordered two-dimensional systems and, 
consequently, conductance distributions of the form given 
by Eq.(fol) are valid in this case. 

Let us now discuss the scaling properties of the mean 
conductance at the band center. Conductance has been 
averaged over a significant number of randomly generated 
samples of increasing linear sizes (Z) keeping almost con- 
stant the total number of measurements. Fig. H shows 
the results obtained by this numerical procedure. The 
scaling of the mean conductance can be fitted to a model 
of the form: 



<G>= 



bU 



(22) 



with a very good precision. The asymptotic value cor- 
responding to the infinite limit is 0.245. Finite values 
of the conductance of two-dimensional systems show- 
ing chi.r.aj. symmetry have been predicted in a number of 
papersQLlEjOlia. Although, the present numerical simu- 
lation nicely supports these theories some caution must 
be used for two main reasons. First, a somewhat practi- 
cal reason. Results heavily depend on the fact that the 
Green function of this problem is just given by only one 
particular state. While the energy of this state is well 
defined theoretically, it could be difficult to make an ex- 
periment at just a particular energy. Previous authors 



on the subject have clearly shown that chirality is lost 
as soon as E = is lcfttl. Second reason is a bit more 
technical. Scaling properties have been obtained for clus- 
ters of an odd number of sites and, therefore, a state at 
E = {). Present computational facilities do not allow to 
prove that scaling of clusters of an even number of states 
proceeds in the same way. (Note that in this case the 
Green function should be recalculated for any new posi- 
tion of the point contacts since storing of the whole Green 
function matrix of the isolated cluster is not possible). 
Nevertheless, I have checked that odd/even differences 
are minimal for small systems. 

The universal conductance distribution at the band 
center of a two-dimensional chiral system shown by the 
thick line in Fig.(Q) is probably the most original result 
in the paper. It corresponds to the L — > oo limit of 
Eq.@ which is obtained for C = 0.0471 in Eq.@ but 
does not differ very much from the distributions corre- 
sponding to the simulated sizes (L = 599, 699, 799) plot- 
ted as three thin indistinguishable lines. All three are 
accurately described by Eq.dTol) with C = 0.0948. Al- 



though obtained for point contacts, previous experience 
with wide wires proves that minor differences should be 
expected for wider but finite contactsllj. Let me insist in 
the idea that perfect agreement with the theoretical con- 
ductance distribution is also explained by a particular 
distribution of wavefunction weight ratios (see Eq.(|20|)) 
that is confirmed by numerical simulation. Fig. @ give the 
distributions obtained for a small cluster size (L = 99) 
-higher maximum value- and for a larger size (L = 499) 
together with its corresponding Gaussian (dashed line). 
The thick line gives the universal prediction for this sta- 
tistical measure of the wavefunction. Consequently, I 
think that unusual phenomena as they have been de- 
scribed in the literature are quite compatjbkj-with present 
results for the conductance distributionQO'LJ. 



VI. DISCUSSION 

The scaling properties of the conductance distribution 
of wide wires showing chiral symmetry has been analyzed 
by numerical solution of a simple model Hamiltonian. 
Predictions of a power law decrease of the mean conduc- 
tance have been confirmed. Besides, the strange scaling 
properties of the distribution give rise to unusual statisti- 
cal properties that have been described in the main text. 
For example, while the slow algebraic decay of the mean 
conductance is due to few but very large fluctuations, 
it has been shown that almost all conductance measure- 
ments are exponentially small for large wires. Actually, 
an exponential law is obtained for the median of the dis- 
tribution and for any generalization of the median based 
on the integral of the distribution. Nevertheless, expo- 
nential decay is not proportional to the sample length 
but to its square root. That is a significant deviation 
from standard localization theory. 

Surprisingly, statistical conductance properties of two- 



dimensional systems connected to point contacts are de- 
scribed at the band center by the same distribution that 
describes wide wires. It has been shown that this statis- 
tics appears due to the existence of an underlying Gaus- 
sian distribution of the logarithm of the ratio of the wave- 
function weights at E = 0. The main difference with 
quasi one-dimensional systems is that now both the mean 
and the median conductance seem to reach a finite limit. 
The universal statistical distributions of conductance and 
wavefunction ratios have been obtained. Although they 
share some properties with the corresponding statistics 
of the Circular Orthogonal Ensemble, it has been shown 



that they differ. This numerical result calls for some the- 
ory able to explain the origin of such a simple statistics 
of the wavefunction at the band center of a chiral system. 
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FIG. 1: Length scaling of the conductance of a disordered 
chain at E — 0. Average conductance is given in units of 
the conductance quantum e 2 /h. Although the variance of 
the conductance distribution is large, the values given for the 
averages show error bars that are sensibly smaller than the 
symbols representing the mean values. The fit at large lengths 
by a power law curve of exponent —0.51 ± 0.01 is also shown. 
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FIG. 2: Length scaling for the model of Fig. hi when the 
conductance distribution is described by an alternative central 
value such as the average value of the logarithm (filled circles), 
the median (empty circles) or the conductance value at which 
the integral of the distribution reaches a value of 0.9 (down 
triangles) . All of them are fitted in the whole range of chain 
lengths by an exponential law decaying as a function of the 
square root of the sample length. 
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FIG. 3: Length scaling of the mean conductance of a wire of 
9 (empty circles) or 19 channels (empty squares) at the band 
center (E — 0). Average conductances are given in units 
of the conductance quantum e 2 /h and fitted to power law 
curves of exponents —0.51 ± 0.01 and —0.54 ± 0.01, respec- 
tively. Filled symbols give the values of the variance obtained 
for the same sets of randomly generated samples. Dashed 
lines are just 2/3 of the previous fits. Error bars are of the 
order of the symbols representing the averages. 
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FIG. 4: Length scaling of the probability distribution of con- 
ductances of a nine channel wire. Filled circles, empty circles 
and diamonds give histograms for the statistics of 10 5 sam- 
ples of lengths 100, 1000 and 10000, respectively. Continuous 
lines give results obtained from Eq. (|l0j). Conductance unit 
is e 2 /h. 
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FIG. 5: Dependence of the mean conductance and variance on 
the number of channels for a fixed length of 2000. The number 
of samples used to get accurate values of the statistical values 
is 10 4 . Two fits to the data are shown: one using Eqs. (4.11a) 
and (4.11b) of the second paper in Ref.ffl with I — 1.646 
(dashed line) and a linear fit (C given by Eq.(n5|)) (continuous 
line). 
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FIG. 6: Length scaling of the conductance of a wire of 9 
(circles) or 19 channels (squares) at the band center (E — 0). 
Scaling of the median (filled symbols) is shown besides the 
data of Fig. H in a semilogarithmic plot. Medians can be 
fitted over the whole range of studied lengths by ([]) giving 
values of £ equal to 5.73 and 8.99 for the 9 and 19 channels 
wires, respectively. 
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FIG. 7: Probability distribution of conductances between 
point contacts in a large disordered square sample. Results 
for three sizes (L = 599, 699, 799) aregiven together with the 
theoretical prediction given by Eq. ([u)|) with C = 0.0948 (All 
four lines are seen as the thin line in the figure) . The univer- 
sal distribution corresponding to the infinite two-dimensional 
system is given by the thick line (C — 0.0471). The conduc- 
tance distribution corresponding to the Porter- Thomas wave- 
function statistics is given by the dashed line. Conductance 
unit is e 2 /h. 
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FIG. 8: Probability distribution of the wavefunction weight 
ratio in a square disordered dot. Results for two sizes (L = 99 
and 499) are given as continuous lines together with the the- 
oretical prediction (Eq. (pJQ|)) for C = 0.09124 which corre- 
sponds to L — 499 (dashed line). The universal distribution 
corresponding to the infinite two-dimensional system is given 
by the thick line (C = 0.0471). 
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FIG. 9: Size scaling of the mean value of the conductance be- 
tween point contacts in a square disordered dot. Continuous 
line is a fit of the form given by Eq.(|22J) with an asymptotic 
value of a = 0.245. Conductance is given in units of the 
conductance quantum e 2 /h. 



